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Abstract and introduction 

Fast Fourier transform methods provide a fast and efficient means of pro- 
cessing large amounts of gravity or geoid data in local gravity field model- 
ling. The FFT methods, however, has a number of theoretical and practic a 
limitations, especially the use of flat-earth approximation, and the requi- 
rements for gridded data. In spite of this the method often yields excellent 
results in practice when compared to other more rigorous (and computational 
lv expensive) methods, such as least-squares collocation. 

The good performance of the FFT methods illustrate that the theoretical 
approximations are offset by the capability of taking into account more data 
irT larger areas, especially important for geoid predictions. For best 
results good data gridding algorithms are essential. In practice trunca 
collocation approaches may be used. For large areas at high latitudes e 
gridding must be done using suitable map projections such as UTM, to avoid 
trivial errors caused by the meridian convergence. . 

In the paper the FFT methods are compared to "ground thruth data in New 
Mexico (5,n from Ag), Scandinavia (N from Ag, the geoid fits to 15 cm over 
2000 km), and areas of the Atlantic (Ag from satellite altimetry using Wie 
ner filtering). In all cases the FFT methods yields results comparable or 
superior to other methods. 

Gravity field modelling by FFT methods 
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The Fourier transformation methods have become practical tools in geodesy 
after the advent of high-degree (180x180) spherical harmonic expansions of 
the geopotential . Using such a geopotential reference model, the residual 
field may with good approximation be treated within the flat-earth approxi 
mation of the Fourier methods. 

For the (residual) anomalous potential T, two-dimensional Fourier trans- 
formation 
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yields the familiar simple frequency domain relationships 
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above equations may be evaluated uy uc rcau — 3 ' . 

the final results obtained by an inverse transform. The use of FFT requires 
gridded data, and introduces errors due to periodicity assumptions and alia 
sing. For details see, e.g., Kearsley et al. (1985). . 

The "inverse Stokes" transformation N to Ag is a high pass filtering ope 
ration, often seriously affected by noise when data is from satellite alti- 
metry'. In this case Wiener filtering may be used. Assuming Kaula s rule to 
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be valid, the PSD of the geoid heights 
1984). In the presence of white noise 
will be 


*NN wili decay like k -4 (Forsberg, 
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The constant c depends on data noise and local variability of the gravity 
field. It may be specified indirectly through the wanted resolution of the 
so u ion, the resolution here defined as the wavelength corresponding to the 
wavenumber with ad^) = 0 . 5 . Resolutions around 20 km give good results. 

T FFT prediction examples in the sequel have been done using the GI 
programme modules GEOGRID (fast gridding using truncated collocation with a 
second order Markov model), GEOPOUR (FFT manipulation ) , and TCIP (interpola- 
tion from result grids with local splines). 


Ex amples of gravity field modelling by FFT 

1. Ag to (5,n) - White Sands area. New Mexico. 

Gravity data was gndded in a 4° x 3° area on a 2* x 2' grid using ter- 
rain reduced data. The FFT prediction results were compared to 384 deflec- 
tion pairs from astrogeodetic observations, yielding 

€" n" 


R.m.s. observed deflections 2.69 6.16 

Difference FFT minus observed 0.73 0.85 


For a comparison r.m.s. prediction accuracies in the range 0. 9-1.0" for £ 

’ 5°L n were ^P 01 ^ in Kearsley et al (1985), using the same 

data with four different prediction methods (ring integration, collocation, 
collocation/integration, and FFT on 5' mean data) 

2. Ag to N for a large area — Scandinavia. 

The geoid was computed for the region 54°-71°N, 4°-32°E by two FET solu- 
tions gridding the available gravity data on either a 6' x 12' geographical 
grid, or a 10 km x 10 km UTM grid (zone 33). The FFT results are compared to 
an earlier comprehensive blocked collocation prediction consisting of 26 3°x 
6 solution blocks (Tscherning and Forsberg, 1986), and all solutions are 
compar to 41 GPS derived geoid undulations (fig. 1) observed by the Insti- 
tut fur Erdmessung, Hannover (Denker, pers.comm. ) . 

_ Th e fFT solutions were both based on GPM2 and a thinned out gravity data 
, b x j" Pixels), as for the collocation solution, with some new data 

a ed in northern Germany and off northern Norway. A total of some 20000 
gravity stations were gridded onto a 210x180 ( geographic ) or 256x256 point 

grid (UTM) . The comparison of the predicted geoids to the GPS results are 
shown below and in fig. 2. 


Geoid comparison at 41 GPS stations (m) : mean std.dev. 


Difference 

FFT (geographic) minus 

GPS 

0.07 

0.17 

— 

FFT (UTM zone 33) - 

GPS 

-0.09 

0.15 

- 

collocation 

GPS 

-1.31 

0.54 

Difference 

FFT (geographic) - FFT 

(UTM) 

0.16 

0.04 


The FFT solutions are clearly superior to the collocation solution, illu- 
strating the benefit of being able to take more data into account. The UTM 
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3. N to Ag - satellite altimetry off Norway and Puerto Rico trench. 

In this example Wiener filtering with resolution at 20-30 km is 
used. With the simple FFT approach useful Ag values were not obtained. In 
both cases the satellite altimetry data was submitted to a local bias cross- 
over adjustment prior to gridding and FPr. 

The test area S of Lofoten, Norway, is a shelf area. Available SEASAT da- 
ta within a 5°x 8° were gridded at 3'x 6' using GPM2 as reference field. The 
adjusted orbits had r.m.s. cross-overs at 8 cm. The computed gravity anoma- 
lies are compared to recent SK gravimetry in a 2.5°x 4 ° area (fig. 3). For 
comparison a collocation solution (GEOCOL) was also done. 

In the Puerto Rico area combined SEASAT/GEDS3 data were compared to all 
ship gravimetry in a 3°x 3° area over the trench. The r.m.s. crossovers af- 
ter adjustment was 31 cm. The FFT solution was done by a 6'x 6' grid in a 
8°x 8° area (without reference field). Results are shown below, together 
with QSU collocation results (Kadir, Knudsen and Rapp, pers.comm) . 

Ag observed FFT difference Collocation diff. 

Area mean std.dev. mean std.dev. mean stddev. 


Norway 12.2 33.8 -1.0 10.9 1.2 12.8 

Puerto Rico -125.8 130.5 -12.7 15.3 -10.5 15.9 


i IT 1 2° 14“ 
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Fig. 3. Gravity ano- 
malies computed by 
FFT from SEASAT al- 
timetry south of Lo- 
foten, Norway. 5 
mgal contour inter- 
val . Comparison 
gravimetry tracks 
shown with crosses. 


The FFT methods work, theoretical approximations are offset by the abi- 
lity to handle more data. For large blocks proper map projections must be 
used. 
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